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m ■ ABSTRACT 



o 



We examine global incompressible axisymmetric perturbations of a differen- 
tially rotating MHD plasma with radial density gradients. It is shown that the 
standard magnetorotational instability, (MRI) criterion drawn from the local dis- 
^ I persion relation is often misleading. If the equilibrium magnetic field is either 

purely axial or purely toroidal, the problem reduces to finding the global radial 
^ ^ eigenvalues of an effective potential. The standard Keplerian profile including the 

I origin is mathematically ill-posed, and thus any solution will depend strongly on 

^ ' the inner boundary. We find a class of unstable modes localized by the form of the 

O ] rotation and density profiles, with reduced dependence on boundary conditions 

(N 
O 
oo 
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^ . Subject headings: MHD, instabilities, accretion disks 



^ ■ 1. Background 



It is often stated t hat the Magnetorotational Instability (MRI) (jChandrasekharl 1 1961 



Balbus fc Hawleylll99ll . hereafter BH91) in accretion disks is a 'local' instability, i.e. normal 
modes are driven unstable by the local value of the rotational flow shear. Implicit in this anal- 
ysis is the assumption that equilibrium rotation and density vary over a much larger spatial 
scale than the mode wavelength. Although it has been shown that short-wavelength hnear 



local MRI mo des can drive global turbulence in the nonlinear regime (IBalbus fc Hawleylll998 



Hawleyi l200ll ). it is worthwhile to study linear instabilities with large radial extent which 
can contribute more direct angular momentum transfer. In order to study these modes, it 
is necessary to use a more complete treatment, including the radial variations of the equi- 
librium profiles. One arrives at a second-order differential equation, which can be solved for 
the structure of global eigenmodes. 
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A major drawback to the study of global modes, aside from the computational com- 
plexity, is their strong dependence on boun dary conditions, which are often unknown for 
astrophysical systems. Previous work (e.g. iDubruUe &: KnoblochI Il993l : lOgilvie fc Pringle 



19961 . hereafter OP96), has included rigid boundary walls to discretize the mode spectru m 



but this imposition is arbitrary and unphysical. Another approach (jCurry et al.l Il994l ) is 
to use pressure constraints to define a boundary, outside of which the f i eld is taken to be 
vacuum. Appropriate matching conditions are then used. iKersale et al.l (120031 ) studied the 
global MRI problem with inflow, and found that certain boundary conditions gave rise to 
"wall modes" with large growth rates. In this paper, we show that unstable axisymmetric 
modes in cylindrical geometry can be described by an effective radial potential. The sign 
of this potential on the boundary dictates how strongly the mode structure depends on the 
specific boundary conditions taken. We find that smoothly varying equilibrium rotation and 
density profiles can localize modes and reduce dependence on the explicit treatment of the 
boundary conditions. The dispersion relation for these global modes can differ greatly from 
that of the local treatment. 

Much previous analytical work on global modes has been carried out using the Boussin- 
seq approximation, which treats both the equilibrium and perturbed density as constant 
except in the equation of motion. This greatly simplifies the dispersion relation, and the 
density gradient appears only through the buoyancy te rm (the Brunt- Vaisala frequency). 
Although this is an appropriate step in the local analysis (IBalbus fc Hawleylll99ll ). when the 
effective radial wavelength is of the order of the system size, we must allow for significant 
density variations over the region in question. In our analysis the density appears in the 
mode equation in two additional ways; by allowing the local Alfven frequency to change 
with radius, and by introducing terms proportional to the first and second derivatives of the 
Alfven frequency. It is well known tha t density gradi ents in ideal MHD can lead to both 
continuous and discrete Alfven spectra (ISedlaceklll971l). In the absence of equilibrium flow, 
these modes lead to damped surface eigenmodes (I Chen fc Hasegawalll974j ). If a free-energy 
source such as differential rotation is present, they can couple to produce unstable modes. 
We examine how astrophysically relevant density profiles can serve to stabilize modes as well 
as help the appropriate imposition of outer boundary condition. 

This paper is organized as follows: In Section [21 we derive the radial global mode 
equation for incompressible axisymmetric perturbations. By allowing for any rotation and 
density profile, our analysis remains relevant to the global MRI in accretion disk s as well as 
other rotatin g systems such as labo ratory experiments (e.g. iGoodman fc Jill2002[ ) and stellar 
core collapse (lAkiyama et al.ll2003l ). In Section [31 we investigate certain illustrative limits of 
the system represented by this equation. Finally, numerical results are presented in Section 
[U showing how certain rotation profiles can lead to direct localization of these modes. 
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2. Basic Equations 

The equation of motion for an MHD plasma is: 



p ( ^ + (v ■ V)v ) = -VP + —(V X B) X B - 

ot J Ho 



Wliere $g is the gravitational potential and P is the scalar pressure. The magnetic field 
evolves according to Maxwell's Equation: 

— = V X (v X B) (2) 

These equations, along with the divergence condition on the magnetic field, admit a rotating 
cylindrical equilibrium of the form Bq = (0, rfigl'"), -B^l?")) and Vq = (0, ri7(r), 0). In this 
equilibrium, we can write ([T]) as the Euler equation: 

p{rQ^r -V<^g) -VPo + —{V xBo) xBo = (3) 

The specification of a density, magnetic field, and rotation profiles will determine the pressure 
up to a constant. We can then define a (local) adiabatic exponent through the identification 

P = Kp^. 

If none of the equilibrium quantities depend on the height z or angle 6, we can take 
Fourier transforms in the axial and azimuthal directions. Neglecting the perturbed grav- 
itational potential 5$^, the equations for the no rmal modes of Lagrangian perturbations, 



(£ = ^(r)e*^'^"^^"^™^ '^*^) to this equilibrium are ( Frieman fc RotenbergI Il960l : IChanmugam 



19791): 

- + 2zp(vo • V)| - ^(0 = , (4) 

Where 

H^) = V(7pV-^ + (^-V)P) + V-(pOV$, 

+ V(Bo-b) + (Bo- V)b + (b- V)Bo (5) 

+ V- (p^(vo- V)vo-pvo(vo- V)0- 

The perturbation of the magnetic field is b = V x x Bq). For incompressible (V ■ ^ = 0) 
perturbations, eq. (5) can be written as two scalar equations: 

(pul-uD-^i^T = [{pu;l-u\){pul-ij\-2prQn' + 2rQBn'B + pN') (6) 



TTl 
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{pcul - col)- — irir) = -2—{puj^n + nBtOA)^r + { ^ + kl] ^Pt (7) 
r ar r \ r"^ J 

Where Um = uj — uiVL^uja = k^Boz + rnQs, i^T is the total perturbed pressure (gas 
plus magnetic), and N"^ = — ^(rfi^ — V$g) is the Brunt- Vaisala frequency. In the present 
treatment, we allow p to vary significantly over the region under consideration; as long as 
we restrict our analysis to incompressible perturbations, the above equations are still valid. 
Through density variation, the local Alfven frequency can change even when the equilibrium 
magnetic field is constant. The density length scale Lp = {d{lnp)/dr)~^ is taken to be much 
larger than the ion Larmor radius (drift waves are not considered). 



2.1. Axisymmetric modes 



For this paper, we restrict our consideration to axisymmetric modes (m=0). Then 
Eqs. (EI) and can be reduced to a single second-order differential equation in the radial 
coordinate: 



d 

dr 



r dr 



F{r) - 2prnn' + 2r^]B^]'5 + pN^ 



F(V) 



= (8) 



Where F{r) = p{r)uj'^ — u}\, and u}\ = klBl^/ pq. This equation describes the standard 
MRI in the limit p' — 0, and the gravitational Interchange Instability in the limit VL ^ 
0, O' 7^ 0. If there is no equilibrium rotation in the system, the mode equati on is the 
cylin drical form of the well known differential equation for surface Alfven waves (jSedlacek 
197ll ). There is a continuum of stable oscillations at each frequency u = ujA{r) = k • 
Bo/ a/ pop{r); each frequency is strongly localized around the characteristic radius where F{r) 
vanishes. These modes overlap spatially and give damping proportional to 1/t. In addition, 
there exists a discrete spectrum of surface modes with position-indepenent frequency. This 
is the phen omenon of damped resonan t absorption, with weak damping for sharp variations 
in density (iHasegawa fc Uberoil[l98^. Finite L armor radius terms couple these modes to 
the Kinetic Alfven Wave (KAW) (jMahajanlll984j ). which will not be addressed in this paper. 



The mode frequency u enters the the differential equation ([8]) only through F and the 
last term in square brackets. If either B^o or Qb vanishes, only u'^ appears. Since all other 
terms are real, the eigenvalues uj must then be purely real or imaginary ( 1 Chandr asekhar 
19601 ). Although the presence of velocity shear makes the evolution equation non- Hermit ian, 
when restricted to the normal mode solutions of purely axial or toroidal fields, we obtain 
a fully Hermitian eigenvalue problem (provided that the equation is well-behaved at the 



- 5 - 



boundaries). This allows for a significant simplification in the search for global modes. 
Before proceeding, we examine the local limit of the above mode equation. 



2.2. Local Dispersion Relation 

In order for modes to be spatially oscillatory, we must have ^" < 0. If the radial 
variation of the equilibrium quantities is small with respect to the scale of the perturbation, 
we can take an expansion ^(r) = 1 + (3r{r — ro)^/2 in equation ([8]). Solving for /5(r = ro), 
and finding which values of ro make /3r < 0, we obtain: 

Pr-K[^^ F + "F J+;^"rF' 

where here all quantities are taken at their local values. In the Boussinesq limit {F' = 
0, 7^ 0), this reproduces the local MRI dispersion relation of BH91, if we identify Pr 
—k"^ and take kr,kz ^ 1/r. Thus the global analysis and the local analysis agree in the 
appropriate limit. However, this equation provides no indication as to which radius should 
be used when applying this criterion, or what to do if it is satisfied in some places and not 
in others. The local instability criterion can be useful in locating the region containing the 
most unstable mode, which for the MRI generally occurs near the point of greatest shear. If 
the mode has radial extent comparable to the equilibrium variation, the full global analysis 
can lead to results quite different from this local criterion. 



2.3. Effective Potential 



When the equilibrium magnetic field is purely axial, eq. ([8]) only admits modes with 



real uj"^. If the Alfven term F(r) 



pto 



uo\ has the same sign for all r in the domain, we 



can make the substitution y = \/ ±rF^r- For purely growing modes, F[r) is negative for all 
r, regardless of the form of the density profile. We arrive at: 

^2 4p2r]2^2~ 



V{r, u)y = 



V{r,uj) = r 1 



2prm' ^ pN' 



Q{r) 



F 

1 F' 1 F" 
^ 2T 



F2 
_ 1F;2 
4F2 



(10) 



The problem becomes one of f inding the zero energy solut ions of the frequency depen- 
dent "effective potential" V{r,uj) (IMahajan fc KrishnanI l2006l ) . If the potential is positive 
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everywhere, the solutions are monotonic, and it is impossible to construct a global solution 
satisfying both boundary conditions. It is therefore necessary that V{r,uj) < in some 
region in r for a global mode to be possible. There are two distinct ways for this to occur: 



1) V{r, uj) is negative all the way up to one or both of the boundaries of the region under 
consideration. This gives rise to boundary-localized 'wall' modes like the ones found in 
OP96. Any change of the boundary condition will drastically affect the mode structure 
and frequency spectrum. The Keplerian flow profile is always of this type for the inner 
boundary, as we shall see below. 

2) V{r, uj) has a minimum which is less than zero at some radius, but is positive elsewhere. 
In this case, the region of oscillation is localized by the potential well, and the mode is 
spatially evanescent outside the well. The boundaries play a reduced role in the mode 
structure, although they can still be important. The local stability criterion may be 
satisfied over significant portions of the disc yet unstable modes can exist which are 
localized by the effective potential well. The most unstable modes are the ones with 
no zero crossings; these modes tend to have a greater radial extent and thus a greater 
chance of carrying radial angular momentum. 



If the equilibrium magnetic field is purely toroidal, the effective potential becomes: 

Since the equilibrium magnetic field is perpendicular to k, the Alfven term F — > pcu^, and 
the only coupling to the magnetic field is through the equilibrium magnetic shear 2rfisf2'^. 
If the modified Rayleigh criterion 

4fi2 + 2rVLn' - A^2 _ 2rVLB9.'jj/p > 0, 

is satisfied, the potential is always positive for purely growing modes, and the system is 
stable to m = perturbations (i.e. there are no global axisymn aetric MRI m odes). In the 



absence of rotation, this is the Tayler "pinch" stability criterion (jTaylerlll973l ). The current- 
free configuration Qb = /3r~^ is always stabilizing. In what follows, we examine only purely 
axial magnetic fields, and defer consideration of toroidal fields to a later paper examining 
non-axisymmetric disturbances. 



3. Global Solutions 



We begin by investigating various limits of the global mode equations analytically. 



- 7- 



3.1. Rigid Rotation 

If r2 = Qq, B = Bzo, and the density is constant, the mode equation reduces to (Fq = 
u"^ — uja/po = const.): 

drrd/^"- ^^^^ 

allowing a family of solutions 



descr ibing shear Alfven waves in a rigidly rotating homogeneous plasma (iHasegawa &: Uberoi 



19821 ). Here /i is to be interpreted as an effective radial wavenumber obeying: 



f^' = kl (^-p^ - 1 j . (13) 

The values of /i^ are determined by matching the solutions of Eq . [T2]to the imposed bound- 
aries. T his results in a (bounciary-d ependant) discrete spectrum of stable eigenmodes when 



/x^ > dPubruUe fc Knoblochlll993h . 



If the dispersion relation of equation (fT3|) gives a negative value for /x^, the solution 
is a linear combination of the modified Bessel functions /i(|/i|r) and Ki{\ii\r). When the 
rotation frequency is constant throughout the entire domain, there can be no global mode 
satisfying both boundaries, as both solutions to eq. f|T2l) are monotonic. If only a portion 
of the domain is subject to rigid rotation (the effective potential is positive in that region 
but negative elsewhere), the modified Bessel functions provide suitable limiting forms. In 
particular, when either the density or the rotation are small for large r, we obtain ^ 
(vacuum solution). We will use this result in Section H] to provide interior and exterior 
boundary matching conditions for modes localized by the form of the equilibrium profiles. 



3.2. An Exactly Solvable Profile 



For the next limit, we investigate a system with differential rotation for which we can 
find exact solutions. Take constant density p = po, let 



"0 



r 



(14) 



and take the magnetic field to be uniform in the z direction. For this profile 1]' < if a > 0. 
We have deliberately chosen this form so that the result for rigid rotation can be obtained by 
letting a ^ 0, /3 ^ 1. A pure power law is obtained when /3 — > 0. As we have stressed above. 
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any rotation profile may be obtained by specifying an appropriate equilibrium pressure. The 
effective potential equation 



4^2 



+ kl 



^0 



1 



(15) 



may be written in the standard Whittaker form (lAbramowitz &: Stegunlll964j ) 



go ^ 3 



y = o, 



(16) 



where 



Po = kl 1- 



-^0 



On a semi-infinite domain r G (0,oo), it allows there are well-behaved solutions ($ is the 
Kummer Function): 



yn 



Qo 



3, 2^r 



(17) 



This solution requires po > 0, which is satisfied for purely growing modes (all modes if 
(3^0). The eigenvalue condition arises from the need for the displacement and its radial 
derivative to be bounded at both r — > and r —>■ oo; the latter demands that the Kummer 
series terminates. This happens when 



3 

2 



Qo 



-n, 



yielding the dispersion relation, 



1 - 



n 



(2n + 3). 



(19) 



If /? — > 0, we can solve equation (fT9!) for the frequency: 



2(2n + 3) 

From eq. (I20l) . we derive the instability criterion: 

aQl > Kvl {2n + 3) 

where Va = oJA/kz- 



k^aVLQ 



(20) 



(21) 
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The radial quantum number n appears explicitly in the spectral relation. The first three 
modes and associated effective potentials are plotted in Figure [H For small fc^, equation fl2Tl) 
provides a more severe constraint on VL' (measured by a) than the local criterion fc^f ^ < 2rVLVL' 
(BH91). We see that as the radial mode number n rises, the instability criterion becomes 
harder to satisfy, and the higher order radial modes are less unstable. Thus the converse 
of eq. (12T|) with n = can be taken as a necessary condition for global stability. In 
Figure [2], we plot the growth rate versus Alfven frequency for the first three radial modes, 
for kz = 1. We see that when uja is small, we have 7 ~ uja/\/^ for all modes. As the 
magnetic field is increased, the lowest order modes remain the most unstable. The growth 
rate for each mode reaches a peak value of 7 = -y/afio/ (4^2^ + 3) at the Alfven frequency 
uja = Ty/aVLokz/ {4:y/2n + 3). The growth rate then decreases as lja increases. There exists 
a critical magnetic field strength above which each mode ceases to be unstable, with the 
lowest order radial mode persisting to the highest field value. 
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2 4 6 8 10 12 14 



Fig. 1. — The profile fl'^ = a/r, p = const., considered in section [3^ admits a discrete spec- 
trum of eigenmodes. For eacli radial quantum number n, there exists an effective potential 
(solid lines), which supports a global mode satisfying the evanescent boundary conditions 
(dashed lines). As n increases, the potential well gets deeper, and the mode extends over a 
greater area. For the unstable branch, higher n corresponds to smaller growth rate. 
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Fig. 2. — Growth rate vs. Alfven frequency for the first three radial eigenmodes of the 
profile = f,p = po,kz = 1. The lowest order radial mode is the most unstable for all 
magnetic field values. As the magnetic field strength is increased, the growth rate of each 
mode first increases and then decreases until the mode is stabilized. 
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In this case, the instabihty criterion, though both quahtatively and quantitatively differ- 
ent from the "local" criterion, can indeed be satisfied for Q' < 0. The eigenvalue problem is 
well posed; well-defined square integrable eigenfunctions are associated with unstable modes. 
The discretization is entirely defined by satisfaction of the outer boundary condition. If we 
had instead imposed hard boundaries at some radii ri and r2, we would arrive at a different 
spectrum by including the second solution to Whittaker's equation, which blows up near 
the origin. This illustrative case shows the perils of a rotation profile which is unbounded 
near the origin (making the inner boundary very important), as well as one which has con- 
stant density and slow fall-off of rotation out to large radius (making the outer boundary 
important). 



3.3. Keplerian Profile 

When the equilibrium pressure and self-gravitation of the plasma are negligible, we 
obtain the Keplerian case, = fig'"^^^^ (^o ~ GM/Rq), whence the Brunt- Vaisala fre- 
quency goes to zero. This is the case mos t often used in the study of thin accretion disks 



flBalbus fc Hawlevlll998l : lFrank et al.ll2002[ ). The effective potential is 



VKep{r,u;) = V{r,u;) = e - + ± + Q(r,u;), 



where we have defined 

and Q{r) is defined as in the previous section. If the density is constant, Q{r) = 0, and 
Amust then be a positive constant if there is to be a spatially oscillatory MRI mode. The 
effective potential Vxep then has no potential well- it takes its minimum negative value on 
the inner boundary (Case 1 described in Section [2l3|) . If the region contains the origin, the 
point r = is an essential singularity {Vkep ~ A/r^ — oo). As the The solutions to the 
differential equation near this singularity, although bounded, have divergent first derivatives. 
The problem of Keplerian rotation including the central point is mathematically ill-defined 
in the constant density case, and no global MRI mode can be supported. 

We now consider the effect of non-constant density on purely growing global modes (7^ = 
0). We assume that for astrophysically relevant cases, p is bounded and that p'(r) < 
(the convective stability criterion is satisfied). The Alfven frequency becomes a function of 
position, and A can become negative if 7^ > 3uj\/ p (recall that uja = k^BQ^/ ^p^po). This 
leads to the conclusion that the maximum growth rate of the local MRI in a Keplerian flow 
profile is '^max = \/^^A{f)- In the global case, however, a mode with a given growth rate may 
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be spatially evanescent in one region but oscillatory in the other- a potential well is created. 
If A remains negative as r — 0, the modes will be well-defined. Since this criterion depends 
on the value of o;^, there will always be some growth rates for which A > sufficiently close 
to the origin. The complete spectrum will still be ill-defined. 

Density variation also introduces terms to the effective potential (eq. [10]) which depend 
on the gradient of the Alfven term F{r) = pu"^ — If the density has a power-law profile 
in some region (p = po^~") a > 0), we find that the part of V that is entirely due to density 
variation is 



2r2(l + Wr'^y \2 

where W = ti}^/(po7^)- The contribution is therefore positive for unstable modes. For 
small r, we have Q ~ 1/r^. This inhomogeneity induced term does not diverge as fast as 
the Keplerian term, and has little effect on the inner boundary. An exponential density 
drop p ~ e~^/^ gives similar results. We conclude that density variation cannot remove the 
essential singularity that arises due to Keplerian rotation. 

The density gradient term Q{r) can be negative if the density profile is locally linear 
and drops near to zero quickly. If p = po(l — A(r — ro)/a) between ro and some r2 = rQ + a, 
then 

. A 1 /l A 1 

Q{r) 



2al + W - ^{r -ro) \r 2al + W - ^{r - tq) 

which is negative for 

{r-ro)>^il + W)-'f 
This has an interesting consequence for models which use the constant density ap proximation 



over much of the range and then assume zero density outside some boundary (e.g. lCurry et al. 



19941 ). While discontinuity matching may be used, a sharp density drop over a small region 
can lead to very unstable wall modes, since smaller values of W are more likely to have 
negative Q{r) over a wider range. Also, the further out the density drop-off occurs, the more 
negative this term will be (if A is small and ro is large). If the rotation is sub-Keplerian, a 
negative radial density gradient implies N"^ < 0, so N'^/F > for unstable modes, and the 
buoyancy terms can help to mediate this effect. 

We see that density gradients in most cases serve to 'shut off' modes that exist for 
constant density. Both the local Alfven variation and the positive Q{r) terms serve to shrink 
the region over which the effective potential remains negative. Eventually for a given k and 
uja, the density profile becomes so steep that the most unstable mode is no longer supported. 
Since lowering the density has the effect of raising the local Alfven frequency, we are in effect 
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raising the 'average' Alfven frequency for the mode. This phenomenon mimics an effective 
raising of the magnetic field, which is known to shut off MRI modes. 



4. Numerical Results 

In this section, we numerically examine how modified rotation and density profiles can 
localize and discretize the unstable radial modes. For given Q{r) and p(r), we select the axial 
wavenumber kz and fiducial Alfven frequency uja, and use a shooting and matching code to 
find the growth rates and radial structures of the unstable modes. At the inner boundary, 
starting conditions for the shooting routine are obtained by assuming that there is some 
radius below which both the density and the rotation can be taken constant, (p = po and 
fl = Qq). The bounded solution in that region is then the modified Bessel function /i(|pok), 
as discussed in Section 13.11 The outer boundary is handled in a similar fashion. Care is 
taken to choose the outer boundary far enough out such that the resulting growth rate and 
mode structure so obtained by does not change appreciably when the boundary is moved. 



4.1. Modified Keplerian Profile 

We do not attempt to model real accretion disc boundary layers near the inner ob 
iect; to do so would requ ire a full treatment of pressure, accretion inflow, etc. (see, e.g. 



Rogava fc Tsiklauril Il993l ). Rather, we wish to examine the simplest rotation profiles rel- 
evant to global MRI. Since the MRI is driven by strong shear, it makes sense to examine 
profiles that are flat for small radii, and which smoothly transition to rapid fall-off for large 
radii. To this end, we take a general form 

= l + ((r"°r.)W ^ " = '■ < (22) 

For the case s = 3/2, ri = 0, this profile approaches Keplerian for r i?o, but tends to a 
constant Qq at the origin. This corresponds to a simple model of an accretion disk where the 
pressure support is only significant below some radius Rq, the inner motion corresponding 
to that of a rigid rotor. The shear is maximized near r ~ 0.342i?05 as opposed to the true 
Keplerian case where the shear remains unbounded as the origin is approached. The density 
is taken to be constant (we shall examine the effects of varying density below). We normalize 
frequencies to the central rotation frequency Qq, and lengths to the fiducial radius Rq. The 
density and magnetic field are normalized such that vaRq^ = for f a (normalized) =1. 

The effective potential for this profile is shown in Figure [3] for three values of the growth 
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rate 7, and for kz — I, va — ■'^5. Modes with more rapid radial oscillation have smaller 
growth rates- this corresponds to a deeper potential well. When the well is very deep, the 
radial wavelength is small, and a local treatment becomes valid, but only near the bottom 
of the well. The most unstable mode has no nodes (n = 0). Although the local criterion for 
instability can be satisfied at some radii for larger growth rates, no eigenmodes exist with 
these larger growth rates which satisfy the evanescent boundary conditions. 
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Fig. 3. — Effective Potential for tlie rotation profile = ^o/{l + r^/^), p = po- For 
kz = 1, va = -25, we plot V{r,u) for three different values of 7 (= —icu). As 7 increases, the 
potential becomes less negative. When 7 = 0.0954, the potential supports the lowest order 
(n=0) discrete radial eigenmode (dashed line). Even though there is a negative potential 
for larger growth rates, no eigenmode exists which satisfies both boundary conditions; 7 = 
0.0954 is the most unstable mode for these and va values and all higher n modes have 
smaller growth rates.. For 7 above ~ 0.15, the potential is no longer negative anywhere. 
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The dependence of the growth rate on the strength of the background magnetic field is 
quahtatively similar to the case considered in Section [3^ and is plotted in FigureHl For small 
values of the magnetic field, we find that the unstable radial modes are very close together 
in growth rate and lie close to the shear Alfven wave dispersion relation 7 = k^VAz- For the 
parameters given, the most unstable mode reaches a maximum at uja = 0.2, 7^=0 = .105. 
At this field strength, the n = 1 mode has 7n=i = .036, implying that after 1 rotation, the 
n = mode dominates by a factor of ~ ~ 20. For stronger background fields, the growth 
rates diminish, and the higher radial order modes are stabilized. Eventually, the magnetic 
field becomes so strong that even the lowest order radial mode is no longer supported- the 
effective potential is not sufficiently deep to support a radial eigenmode. Also plotted is the 
numerically determined critical stability boundary uja cru for the most unstable mode as a 
function of the vertical wavenumber (in the regime ~ Rq^)- Note that in this range, 
the mode can have a significant radial extent thus the modes we are concerned with have 
effective radial wavenumber kr < kz- Since the critical Alfven frequency rises slower than 
linearly with increasing kz, the critical magnetic field required to shut off the instability 
decreases as kz increases. 
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Fig. 4. — a) The growth rate of the first three radial eigenmodes versus magnetic field 
strength for the rotation profile = r2o/(l + r^/^) at fixed kz = 1. We see that as the local 
Alfven frequency increases, we transition from a shear like mode 7 ~ kzVA to a cutoff. The 
maximum growth rate for this wavenumber is '-/max = -105, occurring at kzVA = -2^0. b) 
The critical Alfven frequency (in units of the central rotation frequency) above which no 
unstable global modes are supported as a function of axial wavenumber. Since the slope of 
this plot is less than linear, the critical field strength vacvu decreases with increasing kz- 
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4.2. Density Variation 

We now take for our density profile a form similar to equation (|22|) : 

P^'^'> = 1 I ^° \/ r > Tp ; p = Po, r <rp (23) 

1 + {{r - rp)/a)i 

This profile, like the rotation profile, was chosen to yield constant density below r^, and tend 
to a power law for large r. As noted above, power law density profiles have a stabilizing 
effect on the global modes. The fall off of density for large radius has the effect of smoothly 
transitioning the effective potential to that of a locally stable vacuum magnetic field (as 
p — > 0, ^2 — ^ kz, and the radial solution becomes Ki{kzr)). 

To demonstrate these effects, we examine the maximum growth rates as the transition 
radius is varied, for fixed a = Rq, q = 2. Figure shows the effective potential for 
the most unstable mode as the density transition point Vp is moved inward. When Vp is 
much larger than the radial peak of the constant density mode, there is little effect on the 
mode, as the density is roughly constant over the region where the mode is oscillatory. As 
the pedestal width shrinks, the effective Alfven velocity increases over the region where the 
mode is nonzero, raising the outer edge of the potential well. For fixed central Alfven speed 
va, the frequency of the mode must decrease so that the well remains deep enough to support 
a mode, and the peak moves inward. For rp below ~ 1, there is no longer a possibility of an 
unstable eigenmode. 
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Fig. 5. — For the profile discussed in Section IT2l we plot the effective potential of the most 
unstable mode for three values of the transition radius Vp = {5,3, 1.05}. When is larger 
than the location of maximum shear, the growth rate remains unchanged (7 = .0955 when 
Tp = 5. As the density falloff moves nearer to this point, the n = growth rate decreases 
(7 = .0233 for Tp = 5), and mode becomes more spatially localized (compare to Fig. [3]). The 
mode is eventually cut off for rp ~ 1. 
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5. Conclusion 

When linear perturbations of an inhomogeneous medium have wavelengths compara- 
ble to the equilibrium variation scale, spatial Fourier analysis becomes suspect, and global 
methods are more applicable. We have shown that the study of long radial wavelength in- 
compressible axisymmetric perturbations of a differentially rotating plasma gives rise to a 
effective potential problem with two classes of boundary dependence. In the first class, the 
effective potential is negative up to the boundaries, i.e. the plasma boundaries are locally 
MRI-unstable. The solutions in this situation will always depend strongly on the type of 
boundary conditions imposed. In the case of pure Keplerian rotation, the eigenmode equa- 
tion has an essential singularity at the origin. Physically, this means that the rotational 
shear is maximized on the inner boundary, making the most unstable modes "wall" modes, 
discretized by the imposed boundary. The second situation arises when the equilibrium 
profile is such that local MHD-stability holds at the boundaries. This can happen if the 
rotation shear vanishes for small radius, such as for boundary layer near the central object 
of an accretion disc. The spatial region over which the unstable modes exist are limited by 
the equihbrium flow and density profiles, leading to reduced dependence on the boundaries. 
The depth of the potential well is a decreasing function of growth rate. As the local limit 
is achieved when the potential well is deepest, this result suggests that global modes may 
be more unstable. In addition, growth rates depend on the background magnetic field in 
a complex fashion. For a given perpendicular wavenumber, there exists a critical uja above 
which no unstable linear modes are supported, as in the case of the local MRI, but the cutoff 
values depend on the global properties of the density and shear flow profiles. 

We have not considered the stability of the system to global non-axisymmetric pertur- 
bations, as the effective potential treatment not as readily applicable when m 7^ 0. The 
potential becomes complex, and overstable convective modes can occur. These may have 
larger growth rates than the axisymmetric modes considered above. For similar reasons, 
the case of mixed toroidal and axial fields was not considered. In resistive magnetofiuids, 
complex ax isymmetric disturba nces can manifest as Helical MRI modes which convect along 



the z axis (jRiidiger et a. 



2006h. Ideal MHD He lical MRI for some simple equilibrium pro- 



files were considered in ICurry fc Pudrita (Il995l ). but these modes depend strongly on the 
boundary conditions. An extension to the potential theory described above will be used to 
study the existence and structure of these modes in a forthcoming paper. 



The Authors wish to thank Dr. Richard Hazeltine and Dr. J. Craig Wheeler for useful 
discussions. 
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